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We explore theoretically the optomechanical interaction between a light field and a mechanical 
mode of ultracold fermionic atoms inside a Fabry- Perot cavity. The low-lying phonon mode of the 
fermionic ensemble is a collective density oscillation associated with particle-hole excitations, and 
is mathematically analogous to the momentum side-mode excitations of a bosonic condensate. The 
mechanical motion of the fermionic particle-hole system behaves hence as a "moving mirror." We 
derive an effective system Hamiltonian that has the form of generic optomechanical systems. We 
also discuss the experimental consequences the optomechanical coupling in optical bistability and 
in the noise spectrum of the system. 
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Optomechanics, the study of the mechanical effects of 
light on mesoscopic and macroscopic mechanical oscilla- 
tors, is an emerging field with considerable promise for 
applications ranging from quantum metrology at or be- 
yond the standard quantum limit 0, 0], and from fun- 
damental studies of quantum mechanics, including the 
cooling of macroscopic objects to their quantum mechan- 
ical ground state Q to the study of the quantum-classical 
interface and the generation of quantum superposi- 
tion of macroscopic objects 0-0 • The cornerstone of 
these developments is the control of mechanical degrees 
of freedom by optical forces, most prominently radiation 
pressure. Optomechanical systems of particular interest 
include optical resonators with a movable end- mirror 8j , 
membranes or mirrors placed inside an optical cavity Si- 
ll | , and nanoscale cantilevers [l2| • 

In complementary work, recent experiments have ex- 
ploited the coupling of the center-of-mass degrees of free- 
dom of ultracold atomic ensembles to the optical field in- 
side a cavity to realize optomechanical systems. In such 
experiments an analog of a moving mirror can be pro- 
vide d, e . g., b y side-mode excitations of a bosonic conden- 
sate 3, l3 and by the vibrational motion of cold ther- 
mal atoms [la, Il6|. In a sense, these experiments can be 
thought of as " bottom up" realizations of optomechani- 
cal systems, in contrast to the "top down" systems real- 
ized by mirrors, membranes or cantilevers. Here, photon 
recoil results in the excitation of atomic density fluctua- 
tions associated with the occupation of momentum side 
modes. Experimental signatures of this effect include the 
onset of optical bistability in the transmitted light Q. 

This Letter extends these ideas to ultracold fermionic 
gases, and show how a mechanical mode of ultracold spin- 
less fermionic atoms trapped inside an optical resonator 
can couple optomechanically to the cavity field. In the 
limit of low photon numbers, or equivalcntly the regime 
where only the lowest diffraction order of the atoms by 



the light fields is significantly excited, we identify a com- 
plete analogy between a mechanical fermionic mode con- 
sisting of a superposition of particle-hole excitations and 
the more familiar momentum side mode excitations of 
a bosonic condensate. Evidence of the optomechanical 
coupling can again be observed in the bistable behavior 
in the semiclassical stationary state of the system, as well 
as in its noise spectrum. 

The quantum statistics of atoms are often irrelevant in 
cavity QED experiments using atomic ensembles. How- 
ever, the situation is different when dealing with the 
center-of-mass degrees of freedom of ultracold atoms. In 
this regime of cavity QED and atom optics experiments, 
bosonic condensates are often preferred due to their for- 
mal similarity with photons, and also because one often 
wants to exploit Bosc enhancement. On the other hand, 
condensate dynamics is often dominated by atomic col- 
lisions and collisional frequency shifts that can mask or 
even suppress the effects of atom-field interactions. This 
is one reason for the interest of fermions instead of bosons 
in high-performance atomic clocks [l7, 18 1 and interfer- 
see also Refs. [20-251 which consider some 



ometers [19(, see also Rcfs. [2C 
aspects of fermionic atom optics. 

Figure QJa) shows the system under consideration. A 
Fabry-Perot cavity of length L containing N trapped 
spinlcss fermionic atoms of mass M is driven at rate r\ 
by a pump laser of frequency wl and wave number K. 
We assume that the magnetic trapping frequency is much 
tighter in the transverse y and z directions than along the 
cavity axis x so that low-energy excitations of fermions 
occur predominantly occur along x. We quantize the 
atomic motional degree of freedom along that axis and 
describe its dynamics within a one-dimensional model. 

When the pump laser frequency u>l is far detuned from 
the atomic transition frequency oj a , the excited electronic 
state of the atoms can be adiabatically eliminated and the 
atoms interact dispersively with the cavity field, taken to 
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FIG. 1: (Color online) (a) A sample of two-level fermionic 
atoms with resonant frequency uj a is trapped in a Fabry-Perot 
cavity of length L, and interact with the standing-wave light 
field. The left mirror is partially transmissive and the cavity 
is driven by a laser of frequency lol- (b) Ground-state atomic 
momentum distribution in the absence of grating, with kF the 
Fermi momentum, (c) Atomic momentum distribution result- 
ing from the diffraction of atoms off the grating for low intra- 
cavity photon numbers. States within the width Sk close to kF 
acquire a momentum boost 2K due to photon recoil, with a 
conjugate process close to — kF with momentum boost —2K. 
(d) Energy dispersion corresponding to (c). For 2K <C &f 
the quadratic dispersion relation of the atoms can be approx- 
imated by its gradient at the Fermi energy. 



be single mode. In the dipole and rotating-wave approx- 
imations, the atomic part of Hamiltonian is then [26| 



dx&{x) [p 2 x /(2M) + hU cos 2 (Kx) e'e] (1) 



where i£(x) is the atomic field operator, and c the an- 
nihilation operator of a cavity photon. Because of the 
absence of low-energy s-wave collisions between identical 
fcrmions, they interact only with the light field with the 
coupling Uq = 9q/(wl— w a ), where go is the single-photon 
Rabi frequency. 

We expand the atomic field operator in terms of plane 



waves = i _1/2 Efc/fce 4 

tion operator, with {fk,fli} = 
yields the second-quantized Hamiltonian |22|-|24| 



; where fk is a annihila- 
fa,k', {/fc,AiL=0. This 



k 

+\hU $cJ2(fl + 2Kfk + ftfk+uc) (2) 
k 

where e(k) = h 2 k 2 / (2M) is the kinetic energy of a 
fermion of momentum k, A = uj c — u>l + UqN/2 is the 
effective cavity detuning, and w c the resonant frequency 
of the empty cavity nearest to the laser frequency. 

Because of photon recoil by the off-resonant standing- 
wave field, the atoms can suffer a momentum kick ±2£K, 



I being the diffraction order. For feeble optical fields it is 
sufficient to consider the lowest diffraction order, 1 = 1, 
in which case the physics of the problem is best captured 
by introducing the operators 



(3) 



where p = 2K . It is easily verified that pj, = p- v - 
Equation ([3]) describes a superposition of particle-hole 
excitations with a well-defined excitation momentum p, 
assumed to be positive without loss of generality. We 
further define "right-" and "left-" propagating atomic 
density-fluctuation operators as pp + \ pp \ correspond- 
ing to the summation of k in Eq. ([3]) for k > 0, and k < 0, 
respectively (27j . The examination of the commutation 
relations of p±J leads naturally to the introduction of a 
new set of operators 
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where (3 p = \/2ix jpL = tx / KL is a normalization con- 
stant. These operators obey the bosonic commutation 
relation [6± P ,&± P ] = 1, and [b± p ,b± p ] = [b± p ,b' Tp ] = 0. 
The components with opposite sign thus behave as if de- 
scribing distinguishable oscillators. 

We consider the simple situation where the intracav- 
ity field can be approximated as a plane wave, so that 
only momentum transfer along the cavity axis is consid- 
ered and an effective one-dimensional description is ap- 
propriate. We also assume that the fermions are perfectly 
degenerate and occupy all momentum states with k £ 
[— kp,kp] with interval 2tt/L [Fig. [IJb)], and denote by 
Sk the width of momentum states that can be diffracted 
by the optical field. Clearly that width is at most 2£K , or 
2K for I — 1. The excitation energies of the particle- hole 
excitations at the edges of the interval [kp — 8k, /cp] are 
e(k F -5k + 2K) - e(k F - Sk) = 2h 2 K{K + k F - 8k) /M, 
and e(k F +2K)-e(k F ) = 2h 2 K(K+k F )/M, respectively. 

We assume in the following that the number of atoms 
is large enough that K < k F = irN/L, the opposite of 
the regime k F <C 2K of Ref. [23j]. In that case we can 
approximate the quadratic energy dispersion by its lowest 
order expansion about the Fermi energy, J2k e (^)flfk — 
J2k v F?i\k\flfk, where v F = Hk F /M is the Fermi velocity, 
see Fig.QJd). This amounts to neglecting the dependence 
on k of the excitation energies of the particle-hole pairs, 
K(k F ± K) ~ Kk F . From the commutation relations of 
b operators and ^ fc Hv F \k\fl fk, one can then rewrite as 



J2 e ( fc )/fe A -> J2 hv FP( b l b P 



V- p b-p) 



(5) 



p>0 



and the final form of the effective Hamiltonian is 
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H eS = hw M (bl K b 2 K + V_ 2K b-2K) + h A + g{b\ K + b 2 K) + g(b- 2 K + b-2x) c T c + ihn(c ! - c) 



(6) 



This Hamiltonian bears a marked similarity to the 
generic Hamiltonian of cavity optomcchanics, with me- 
chanical oscillator frequency lom = 2Kvf, and optome- 
chanical coupling g = Uq/(AP 2 k)- It is also formally 
identical to the effective Hamiltonian that describes op- 
tomechanics experiments using a bosonic condensate [13J , 
except that now it is the collective density fluctuation 
of the fermions that plays a role of a moving mirror. 
In contrast to the situation with a Bose condensate, in 
the fermionic system the noninteracting ground state is 
a filled Fermi sea. The lowest momentum side mode cor- 
responds then to the process |A;| — > \k + 2K\ rather than 
— ^ |2A'|. A key observation, though, is that although 
the original statistics of the atoms is fermionic, the exci- 
tation is bosonic, due to the fact that a single photon is 
associated with each specific particle-hole pair. 

We next show that there exists an experimentally 
achievable set of parameters that satisfies the conditions 
invoked in deriving the effective Hamiltonian. To guide 
this discussion we concentrate first on the bistable be- 
havior in the system's steady state, one of the hallmarks 
of optomechanical coupling. We proceed by introduc- 
ing the quadratures of the mechanical oscillators Xm 



, P M = P++P- where X± = {b\ 2K +b 

P± = i(b±2K ~b±2K)/V%, respectively. The correspond- 
ing Hcisenberg-Langevin equations are obtained from 
Eq. © as 



dXj£ 

dt 
dc 

dt ~ 



umPm, 



dP M 
dt 



-OJM 



X M - 2V2gc i c, 



A + V2gX 



M 



C + T) — KC+ v2/tCi n , (7) 



where k is the cavit y de cay rate and Ci n denotes a Marko- 
vian noise operator (28| associated with a reservoir at zero 
temperature responsible for the cavity decay. It has zero 
mean and correlations (c\ n (t)c m (t')) = 2n5(t — t'), and 
(ci n (t)ci n (t')) = 0. In contrast with the situation in op- 
tomcchanics with real mirrors, thermal damping can be 
ignored here. 

The mean-field steady-state solution of the 

Langcvin equations (J7J is readily obtained by setting all 



time derivatives be zero. 
-2^2g\c^\ 2 /LU Ml and 

= 
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ft + z (A — 4g 2 ui M \c 



(8) 



This form of cubic equation is characteristic of optical 
multistability 0, [IH [l5 , 25 1 . We consider for illustration 
the experimentally achievable parameters Uq = 2~k x 20 



kHz, A = 500 nm (K ~ 10 7 m" 1 ), L = 100 /xm, N ~ 5000 



so 

25 



atoms yielding a Fermi momentum kp ~ 10 m , 
that k F = 12.5K, k = 2ir x 1 MHz, M = 1.5 x 10 
kg, and the Fermi frequency ujf = £F/h~ — 10 MHz. 
This value of Uq assumes a single-photon Rabi fre- 
quency go = 2tt x 10 MHz and the pump-atom detuning 
bJL — oj a = 2tt x 30 GHz, similarly to Ref [13j. Since the 
decay rate of the atomic excited state is of the order of 
MHz, spontaneous emission can safely be neglected here. 
Typical frequencies of the side-mode excitations are then 
h(Sk) 2 /(2M) ~ h{2K) 2 /(2M) ~ 0.1 MHz, so that the as- 
sumptions invoked in deriving the effective Hamiltonian 
are satisfied. 

Figure [5] shows the steady-state mean photon num- 
ber as a function of the pump-cavity detuning (a) for 
fixed values of tj/k, and (b) as a function of the pump 
rate for a fixed value of the pump-cavity detuning. From 
the linear stability analysis described below, the steady 
state is easily seen to be bistable. The fluctuations 
around the steady state can be studied via a lineariza- 
tion of the Langevin equations about their steady-state 
mean. We proceed by expanding each operator around 
its steady-state value as 0(t) = + SO(t), and intro- 
duce the cavity-field quadratures SX = (Sc* + 6c) /y/2, 
8P = i(5c> — 8c)/V2. The linearized Langevin equations 
are then written in the usual form f(t) = Jf(t) + 
where f(t) = [8Xm,SPm,SX,8P} t , the noise term 
£(t) = [0,0, V2ndX m , V2n8P m ] T , and the drift matrix 
J is given by 



L0 M 

-ujm -4#c s 

-K A 

-2gc s —A -K 



(9) 



Here we have taken c^ s ' to be a real number without 
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FIG. 2: (a) Steady-state intracavity photon number as a func- 
tion of (a) pump-cavity detuning for - by increasing maxima 
- tj/k = 1.2, 2.4, 3.6 (bistability threshold), 4.8, and 6.0; (b) 
pump rate for ljl — w c = 2tt x 47.5 MHz. 
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2 is the effective 



loss of generality, and A = A 
detuning. 

The eigenvalues of the matrix © are obtained by di- 
rect diagonalization, and we find that the steady state is 
stable only if the real part of the eigenvalues is zero or 
negative, and as expected, the branch of the bistability 
curve with negative slope (dashed part in Fig. 2) is found 
to be unstable. Further, the imaginary part of eigenvalue 
corresponds to the oscillation frequency of each quadra- 
tures. At the crest of the bistability curve, the imaginary 
part becomes zero. 

The atom-cavity system is typically probed via mea- 
surements on the optical field transmitted by the cavity, 
in particular its spectrum. The noise spectrum of an ar- 
bitrary stationary field with fluctuations Sf[t] is defined 
by the Fourier transform of its autocorrelation function 
(6f[u]6f[ui']) = S f [w]6{u + u') as 

5/M = h {/ ^'e- i(wW)t (£/>]£/[u/] 

where 6f[us] is the transform of 6f[t] and {■■■}* denotes 
a time average. With 5X in = (Sc- m +6c- m ) / y/2 and 5P\ n = 
i{Sc\ n — <5cin)/\/2, the solution of the linearized Langevin 
equations gives 



6X M [u>] 



4V2k0c< s W{(k - iuj)5X in [Lo} + A5P in [u]} 



d[u 



(10) 



where d[ui] = (uj 2 -u} 2 M )[(K-iu>) 2 + A 2 ]+2uj M A(2gc^) 2 , 
so that 



S 



2/ t (4< ?C ( s W) 2 



(A + lu) 



MM 



(ii) 



and the measurable quadratures of the optical field [2£ 
are 



,5' 



(2gc s YA 2 Sx M + 2k 



-V, 



■A) s 



■5/ 



MM I 2 

(2gc s ) 2 (n 2 + ui 2 )S Xm + 2k \k 2 + (uj + A) 



(12) 



\d[u 



In conclusion, we have discussed theoretically the pos- 
sibility of realizing optomechanical systems using a de- 
generate gas of spinless fermions. We found that a me- 
chanical fermionic mode analogous to the conventional 
moving mirror of optomechanics is provided by the col- 
lective density fluctuations associated with particle-hole 
excitations, and this mode can be quantized as a bosonic 
excitation. We have also identified experimentally real- 
istic parameters such that this effective fermionic mirror 
results in a bistable steady-state behavior of the intracav- 
ity photon number. Generally speaking, optomechanics 



with fermionic mirrors should be observable under con- 
ditions similar to the experiments using bosonic conden- 



sates 



Future studies will include a self-consistent determi- 
nation of the fermionic state in quantized light fields, 
transverse effects, higher-order atomic diffraction, as well 
as the role of the trapping potential, which modifies the 
mean atom number at each site of the optical potential 
and may lead to parametric heating. 
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